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Abstract: Indirect detection signals from dark matter annihilation are studied in the positron channel. We 
discuss in detail the positron propagation inside the galactic medium: we present novel solutions of the 
diffusion and propagation equations and we focus on the determination of the astrophysical uncertainties 
which affect the positron dark matter signal. We show that, especially in the low energy tail of the positron 
spectra at Earth, the uncertainty is sizeable and we quantify the effect. Comparison of our predictions with 
current available and foreseen experimental data are derived. 



Introduction 

The quest for the identification of dark matter 
(DM), together with the comprehension of the na- 
ture of dark energy, is one of the most important 
and potentially far-reaching problems in the un- 
derstanding of the physical world, and represents 
a matter of investigation where cosmology, astro- 
physics and particle physics converge together. It 
is therefore of utmost importance to address the 
problem of the identification of DM from many 
different points of view, which bring together re- 
search performed, from one side, in particle accel- 
erators and on the other side in underground labo- 
ratories, in neutrino telescopes, in large-area sur- 
face detectors as well as in space. In this paper 
we attack the problem from the point of view of 
the indirect detection of DM through the search of 
its products of annihilation inside the galactic halo. 
We will specify our discussion on the positron sig- 
nal, for which we intend to quantify the astrophys- 
ical uncertainties and derive the ensuing conse- 
quences for upcoming experiments. A more ex- 
tended and detailed discussion of our analyses will 
be presented elsewhere [3]. 



The diffusion equation and its solution 

The propagation of positrons in the galactic 
medium is governed by the transport equation: 

-K(E)V 2 *(E)--^{b(EME)} = Q(x,25), 

where ^>(E) is the positron number density per 
unit energy and Q(x, E) is the positron source 
term. The transport through the magnetic turbu- 
lences is described by the space independent diffu- 
sion coefficient K (x, E)—K a e s where e — E/E 
and Eq = 1 GeV. Positrons lose energy through 
synchrotron radiation and inverse Compton scat- 
tering on the cosmic microwave background radi- 
ation and on the galactic starlight at a rate b(E) = 
Eq c 2 /te where te = 10 16 s. The source term 
for positrons produced through DM annihiliation 
is then: 



Q(x,E) 



(av) 



f(E) p 2 (x) =«/(£) 



P 2 (x) 
Pi ' 



where m x is the mass of the DM particle and (av) 
is the low-energy velocity-averaged DM annihila- 
tion cross-section. The positron spectrum f(E) 
produced in the annihilation process depends on 
the specific final states of the annihilation cross 
section. The DM density profile of the galactic 
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Figure 1 : The halo function X is plotted as a func- 
tion of the positron injection energy Es, for an ob- 
served energy of 10 GeV. The galactic halo pro- 
file is NFW [9] with a scale radius of 20 kpc. The 
curve labeled as MED corresponds to the choice of 
propagation parameters which best-fit the B/C ra- 
tio [7]. The MAX and MIN configurations corre- 
spond to the cases which were identified to produce 
the maximal and minimal DM antiproton fluxes 
[4], while the larger area denotes the full variation 
of the propagation parameters yet compatible with 
the B/C analysis [7]. 

halo is denoted by p(x) and its solar value by p Q 
whereas a = 1/2 (Majorana) or 1/4 (Dirac). 

The transport of positrons through the Galaxy 
takes place into a diffusion zone (DZ), which we 
model as a cylinder with medium height L z and 
radius R g = 20 kpc. Until now, the diffusion equa- 
tion has been solved using the method of Green 
functions and we have improved that technique 
with the help of Bessel functions which allow 
to take naturally into account the galactic radial 
boundaries (for details, see [3]). The positron den- 
sity may be expressed in terms of a halo function 
I as: 

roo 

*(E) = k / dE s f(E s ) 1(Xd). 

Eoe 2 Je 

The typical positron diffusion scale Ad depends on 
the injection energy Es and on the observed en- 
ergy E as shown in [5]. That halo function may be 



written as: 

X{\ D ) = Y. C 9 X 3 (x )exp (-^Xl) , 
where 

C ' JUL* * P% X9(X) ' 

and Xg are eigenfunctions of the Helmholtz equa- 
tion [3], The halo function X is featured in Fig.l 
as a function of the injection energy E$ for a fixed 
value of the observed positron energy E and as- 
suming a NFW density profile [9]. In order to 
quantify the uncertainty related to the properties of 
the diffusion region, we have varied the astrophys- 
ical parameters which govern the diffusion process 
inside the range allowed by the B/C measurements 
[7]. Fig.l shows that, for large injection energies as 
compared to the observed energy, the uncertainty 
on the halo function is sizeable, reaching a factor 
of two upwards and of 20 downwards at Es = 1 
TeV. Other DM profiles will be discussed in [3]. 

The positron signal 

The solution of the diffusion equation allows us 
now to calculate the positron energy spectra at the 
Earth, and to determine the uncertainties on the 
fluxes. One example is given in Fig. 2, where the 
positron flux for a DM particle which annihilates 
purely into a fr final state is shown. The annihi- 
lation cross section, here and in the following, is 
chosen to match the relic abundance for cold DM, 
as determined by the latest WMAP data [10], un- 
der the hypothesis of dominant s-wave annihila- 
tion for the DM candidate. The figure shows the 
band of allowed values for the positron signal in 
this channel: the uncertainty ranges from a factor 
of 5 at energies around 1 GeV down to a factor of 
2 above 10 GeV. 

A comparison with the HEAT data [2], in terms of 
the positron fraction, is given in Fig. 3, for anni- 
hilation into a W + W~ final state by a 100 GeV 
DM particle. Due to the contribution of the back- 
ground [8] (which we have reduced by 10%), the 
uncertainty on the predictions for the total flux is 
limited to a factor of the order of 20-30%. Clearly 
a complete understanding of the fluctuations of the 
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Figure 2: Positron fluxes for a 100 GeV DM par- 
ticle which annihilates into a fr final state, as a 
function of the positron energy. The galactic halo 
profile is NFW [9] with a scale radius of 20 kpc. 
The upper, median and lower curves show the ef- 
fect due to galatic propagation: the median curve 
refers to the best-fit set of propagation parameters 
[7], while the upper and lower curves delimit the 
uncertainty band on the theoretical determination 
of the flux. 

theoretical predictions on the total positron frac- 
tion arising from astrophysical origin will require 
a study of the background in the same theoretical 
framework adopted for the signal. This would al- 
low also to correlate the determination of the two 
positron components in the comic rays. This is be- 
yond the scope of this paper, which concentrates on 
the determination of the uncertainty arising specif- 
ically from the signal, and will be addressed else- 
where. From the same figure we also notice, as 
shown by other authors, that the HEAT excess 
around a positron energy of 10-30 GeV could be 
explained by DM annihilations into W + W~ pairs, 
provided that a boosted DM density is present. In 
this analysis, we furthermore reach the robust con- 
clusion that this agreement holds whatever the cos- 
mic ray propagation model. On the other hand, the 
HEAT data alone cannot be used to discriminate 
between the various possible propagation models. 
However, this discrimination will become possible 
with the upcoming PAMELA [6] experiment. The 
predictions for a 3 years flight, in case of a HEAT 
excess explained by a 100 GeV DM annihilation 
line, are presented in Fig.4. Should 
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Figure 3: Positron fraction e + /(e + + e~) as a 
function of the positron energy. The lower solid 
curve refers to the background estimate of Ref. [8], 
rescaled down by a factor 0.9. The upper, median 
and lower curves show the effect of including a sig- 
nal component calculated for a 100 GeV DM parti- 
cle which annihilates into a W + W~ final state and 
for a NFW halo profile [9] with a scale radius of 
20 kpc. A boost factor of 32 is included. The me- 
dian curve refers to the best-fit set of propagation 
parameters [7], while the upper and lower curves 
delimit the uncertainty band on the theoretical de- 
termination of the flux. The experimental points 
refer to the combined 1994, 1995 and 2000 HEAT 
results [2]. 

PAMELA detect a feature like the one shown in 
Fig.4, this would be a clear signature of a DM sig- 
nal (a direct production of positrons, in this case) 
and moreover the experiment would have poten- 
tially the capability to distinguish among differ- 
ent propagation models. Similar conclusions, with 
smaller error bars, will be obtained by the AMS 
experiment [1] on the International Space Station. 

Conclusions 

In this paper we have discussed the positron in- 
direct detection signal of DM annihilation in the 
galactic halo, focussing our attention to the deter- 
mination of the astrophysical uncertainties on the 
positron flux due to the positron propagation in- 
side the galactic medium. We have shown that the 
uncertainties may be large: for a 100 GeV DM par- 
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Figure 4: Predictions for the positron fraction 
e + /(e + + e - ) as a function of the positron en- 
ergy for the PAMELA flight and for a DM panicles 
which annihilates directly into e + e~. The lower 
solid curve refers to background only [8] . The me- 
dian curve refers to the best-fit prediction, while 
the upper and lower curves show the results for the 
extreme values of the propagation parameters. The 
error bars denote the predictions for a detector ac- 
ceptance of 20.5 cm 2 sr [6] and a 3 years flight. 

ticle, they are of the order of 10-30 for a few GeV 
positron, and a factor of 2-5 above 10 GeV, de- 
pending on the propagation parameters, the halo 
profile and the injection spectra. The compari- 
son with current data shows that agreement be- 
tween the predictions and the possible HEAT ex- 
cess is possible, for DM annihilating mostly into 
gauge bosons or directly into a positron-electron 
pair, and the agreement is not limited by the as- 
trophysical uncertanties. We then showed that up- 
coming experiments, like the PAMELA mission, 
will have the capabilities, in the case the excess is 
confirmed and more accurately determined, to shed 
light also on the astrophysical models underlying 
galactic propagation. 
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